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Abstract 

We present a theory for surface polaritons on ferroelectric-antiferromagnetic materials with 
canted spin structure. Canting is assumed to be due to a Dzyaloshinkii-Moriya interaction with 
the electric polarisation and weak ferromagnetism directed in the plane parallel to the surface. 
Surface and bulk modes for a semi-infinite film are calculated for the case of transverse magnetic 
polarisation. Example results are presented using parameters appropriate for BaMnF4. We find 
that the magnetoelectric interaction gives rise to "leaky" surface modes, i.e. pseudosurface waves 
that exist in the pass band, and that dissipate energy into the bulk of material. We show that these 
psuedosurface mode frequencies and properties can be modified by temperature, and application 
of external electric or magnetic fields. 

PACS numbers: 71.36+c;78.20.Jq;78.20.Ls 
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I. INTRODUCTION 



Magnetic polaritons are electromagnetic waves that travel in a material with dispersion 
and properties modified through coupling to magnetic excitations 1 . Polaritons can display 
a number of interesting and useful properties, including localization to surfaces and edges, 
and non-reciprocity- - —, whereby propagation frequency may not symmetric under direction 
reversal: i.e. u(k) ^ u(—k). Surface polaritons at optical frequencies have received much at- 
tention in recent years, and appear in a number of different applications including detectors^, 
biosensors^ and microscopy^. 

Theoretical treatments for ferromagnetic polaritons^ and simple antiferromagnets^ 1 ^ 
were made several years ago. A most interesting class of polaritons are in multiferroic 
materials where magnetoelectric interactions couple magnetic and electric responses^^2. A 
focus of theoretical work has been on bulk modes in linear magnetoelectric coupled media^ii. 
Surface modes have also been discussed for the case of no applied external magnetic or 
electric fields, and neglecting canting of magnetic sublattices^ 1 ^. 

Modification and control of multiferroic surface polaritons through external electric and 
magnetic fields is an intriguing prospect. In the present paper we discuss in detail how 
temperature, electric and magnetic fields affect surface modes in canted spin multiferroics 
with linear magnetoelectric coupling. We allow for canting of magnetic sublattices, and show 
that canting is very important for understanding and manipulating surface mode frequencies. 
Most significantly, we show that the transverse magnetic field polarisation (TM) surface 
polariton excitations are in fact pseudo-surface modes characterized by complex propagation 
wavevectors whose imaginary parts are proportional to the strength of the magnetoelectric 
interaction. 

Linear magnetoelectric coupling is believed to operate in multiferroics BaMnF^ and 
FeTiO^. In order to make contact with previous work on bulk polaritons, we concentrate 
here on the BaMnF4 system. The paper is organized as follows. The geometry and energy 
density of the system are considered in Section II where we also discuss the canting angle 
in relation to the magnetoelectric coupling. In section III, susceptibilities are derived using 
Bloch and Landau-Khalatnikov equations. The electromagnetic problem is solved in Section 
IV and results given in Section V for surface and bulk modes on BaMn4. In Section VI, 
effects on surface mode properties due to possible modifications of material parameters are 



2 



discussed. Conclusions are given in Section VII. 



II. GEOMETRY 




FIG. 1. Geometry. Canting of two magnetic sublattices {m a and nib) by an angle 9 produces a 
weak ferromagnetism (M) along the x axis parallel to the surface. The spontaneous polarization 
(P) is assumed to lie in a plane parallel to the surface. Propagation of the surface mode is along 
the y axis with wavenumber k y . 



The geometry is sketched in FigJH We consider a semi-infinite multiferroic film that fills 
the half space z < 0. The magnetic component of the multiferroic is a two sub-lattice anti- 
ferromagnet with uniaxial magnetic anisotropy. The two magnetic subdattices are allowed 
to cant in the x — z plane with canting angle, 9. We assume symmetric canting such that 
\m a \ = |m&| = M s . The canting generates a weak ferromagnetism which is perpendicular to 
the spontaneous polarisation. This configuration represents a Dzyaloshinkii-Moriya canting 
driven by spontaneous polarisation. Both the weak ferromagnetic moment and spontaneous 
polarisation are constrained to lie in x — y plane, parallel to the surface. The magnetic 
easy axis is out-of-plane, along the z direction. An external electric field is applied parallel 
to the spontaneous polarisation, and an external magnetic field is applied along the weak 
ferromagnet moment. 

For the polariton propagation, we consider in this paper transverse magnetic (TM) po- 
larization in which the magnetic part of the electromagnetic wave propagates parallel to the 
surface. We consider only surface modes traveling along the y direction, so that the magnetic 
component lies in x direction while the electric component has E y and E z components. 
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A fourth order Landau-Ginzburg energy density is assumed to describe the dielectric 
contribution to the energy: 

Fe = \CiP; + + (P* + P 2 ) + \ A 2 (P* + Pt) - P y E y . (1) 

The first and second terms on the right hand of Eq.flTJ represent the energy density for 
the y component of the polarization with (1 and (2 dielectric stiffnesses. The third and 
fourth terms represent the contribution of x and z polarization components with dielectric 
stiffnesses Ai and A2. The last term is the external electric field applied parallel to the 
spontaneous polarization P y . 

The magnetic contribution to the energy density is assumed to be of the form: 

h \M a ■ zf + (M 6 • z) 2 ] - (M a + M b ) x H . (2) 



? M = \M a ■ M b — — 



The first term on the right of Eq. ([2^) is an exchange energy with a strength A > 0. The 
second term represents the anisotropy energy with anisotropy constant K and the last term 
is the Zeeman energy from an external magnetic field. Although the easy axis is out of 
plane, we assume that the symmetry of the canted antiferromagnet results in a negligible 
magnetization along the z direction and therefore ignore demagnetization effects. 

The weak ferromagnetic moment is denoted M x , and the longitudinal component of the 
magnetization is L z . These are defined by: 

M x = (m a + m b ) x = 2M S sin 9 (3) 

and 

L z = (rh a - m b ) z = 2M S cos 9. (4) 
A linear magneto-electric coupling is assumed of the form 

Fme= -aP y M x L z 

= -aP y y ■ (m a x m b ) - aP y [{m a ) x (m a ) z - (m b ) x (m b ) z ) (5) 

where a is the magneto-electric coupling constant. This energy governs the canting of the 
magnetic sub-lattices. 

The canting angle 9 is determined by minimizing the magnetic and magnetoelectric en- 
ergies with respect to 9. Minimizing, we arrive at the condition 

H cos 9 - \kM s sin 29 + 2aP y M s cos 29 + \M S sin 29 = 0. (6) 



In the absence of an external magnetic field, Eq.([6]) simplifies to 

AaP 

1W20) = \ ■ (7) 
y ! K-2X K J 

Note that a positive magnetoelectric constant describes a weak ferromagnetism M x aligned 
along —x. 

The canting angle depends on the equilibrium magnitude of P y , and this is found by 

dK 



minimizing the dielectric and magnetoelectric energies. Requiring -^-{F E + F ME ) = 



results in 

CiPy + C2Py ~ 1olM 2 s sin 26-E y = 0. (8) 

Lastly, the magnitude of Mg depends on temperature. In mean field, the magnitude can 
be written in terms of the Brillouin function B (77) as 

M s = M s (0)B s ( V ) (9) 

where 

= 9pBS_ r_ AM cog 29 + KMg cog2 g + 2ap gin 2Q + ^ gin 0i / 10 x 
k B T 1 y j 

The spontaneous polarization, magnetization and canting angle are calculated by solving 

simultaneously Eqs.([6]), ([8]) and fl9]). Solution for general angles is done numerically using 

root finding techniques for coupled transendental equations. 

III. DYNAMIC SUSCEPTIBILITY 

In order to solve the electromagnetic boundary value problem for the surface and bulk 
polariton modes, we need constituitive relations for the dielectric and magnetic responses. 
We consider linear response and calculate the permitivity and permeability using equations of 
motion derived from Eqs.([T]),02])and ([5]). The equations of motion for the magnetic response 
are given by Bloch equations, 

M = 7 Mx (^i(F M + F ME )) (11) 

where 7 is the gyromagnetic ratio. The equations of motion for the polarization response 
are given by Landau-Khalatnikov equations 

P=-f®(F E + F ME ) (12) 
oP 



where / is the inverse of phonon mass. 

The set of dynamic equations appropriate for the polarizations of TM modes are, 



— ium x = {uja cos 9 + 2u me sin 9) l y (13) 

— iujly = (2u ex cos 9 — u a cos 9 — 4w me sin 9) m x 

+ (2u ex sin 9 + 4:U me cos 9 — u a sin 9 + uj a ) l z 

+A^aM 2 s cos 29p y + 2^M s h x cos 9 (14) 



uj 2 
/ 



— iul z = — (2u ex sin 9 + 2o; me cos 6 1 + oj ) l y (15) 
■Py = (Ci + 3P 2 C2) P y - 2aM s (m x cos 6 1 + l z sin 6 1 ) - (16) 



and 



j Pz = (A 1 + A 2 P 2 )p z -e z . (17) 

The notation used above is in units of frequency and defined as: u a = ^KM S is the 
magnetic anisotropy, u ex = j\M s is the exchange, u me = jaP M s is the magnetoelectric 
coupling and u = ^H Q is the external magnetic field. 

Equations (fl3"|) to (|T5|) are coupled to (flB]) through ME susceptibilities x me an d X em 
defined as 

m = X m h + x me e and p = x e e + x em /*- (18) 
The relevant magnetic susceptibilities for TM modes are given by 

X™ = K cos 2 9 + u me sin 29) { + ° ex ) (19) 

where u s = ^AttM s . The electric susceptibilities are 



V'« M 2 \ + r ? m 2 J <*)) 

(u 2 ey - u 2 ) (< - w 2 ) 



* = R^ 5 )' (21) 



The magnetoelectric susceptibility is 



me _ em _ , )^rne ^2) 
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The frequencies u ey and u m are denned as co 2 y = u 2 y + 5 and Cj 2 ^ = — 5, where 5 is 



expressed in the form 



\ { [K - « - O } (23) 



,= 2 
where 

= C c u; s cos 26 1 (2a; ex sin 2 9 — u a cos 2 # + w sin 6>) (24) 

with C c = ~a 2 M 2 s f. 

u ey is the frequency of the soft phonon along the spontaneous polarization and u m is the 
magnetic resonance frequency: 

W m — &lfm + + ^o- (25) 

Here u a f m is the antiferromagnet resonance, 

tilfm = u a (u a - 2u ex ) cos 2 9 + 2u ex (2u ex - u a ) sin 2 9 (26) 
Q me is related to the magneto-electric interaction, 

= 8w L + 2w me (w« + 2w ea; ) sin 20 (27) 
and f2 Q is related to the external magnetic field, 

Q 2 Q = u [uj + 6u me cos 9 + (4w ex - u a ) sin 0] . (28) 

The frequencies u a = jKM s , u ex = 7 AM,, u me = ^aP Q M s and u = ^H Q represent 
contributions from energies associated with the magnetic anisotropy, exchange, magneto- 
electric coupling and external field respectively. The frequency u ez is the phonon frequency 
along the z direction. 

Other parameters in Eq. fjTTJ]) to ( I2"2"j) are defined as 

C mx =(u%-u 2 m )/(u%-u 2 m ) (29) 
C ex = 5/{u%-u 2 m ) (30) 

Cey=(u%-Ul)/(U%-U 2 m ) (31) 

C my = 5/(u 2 ey -u 2 m ) (32) 
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and 

C me = C af u s cos 26 (u a cos 6 + 2u me sin 6) I {bj\ y - Ujfy , (33) 
where C a f = -aM s f. 

From the expression of the susceptibilities above, it can be seen that the applied magnetic 
field directly influences the susceptibilities through Q a . By way of contrast, the applied elec- 
tric field changes the susceptibilities indirectly by affecting the magnitude of the spontaneous 
polarisation. 



IV. THEORY FOR BULK BANDS AND SURFACE MODES 

Dispersion relations for the bulk modes are obtained by solving the electromagnetic 
Maxwell equations 

Vx^-if (34, 
c ot 

where the fields D with E and B with H in TM modes are connected through constitutive 
equations 

B x = fi x H x + 47i X me E y , D y = e y E y + 4n X me H x and D z = e z E z (35) 

with permeability and dielectric functions defined as fi x = 1 + 47rx™ and e y = 1 + ^x y - 

Plane waves are assumed for bulk traveling modes of the form E,H ~ e i(k y y+k z z-u>t) _ 
Substitution into the Maxwell equations provides an equation for the bulk mode dispersion: 
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The dispersion relation has solutions determined by the zeroes of the dielectric constant 
e z , zeroes of the function f(fi,e) = fi x e y — (47rx me ) 2 . Equation (136)) diverges at the pole of 
dielectric constant e z , e y and x me - 

The dispersion relation for surface modes is calculated by assuming surface localized plane 
wave solutions of the form: 

E, H ~ e^e^yv-^) for z < (37) 



S 



and 

E, H ~ e -Poz e i{kyV-**) for 2 > 0. (38) 

Substitution of Eqs. f|3T|) and (|38p into Eqs. (13~4"|) provide an implicit relation for the 
attenuation constant in the material, 

e z [p + %^ X me -^) = e y k 2 y - e y e zf x x {^)\ (39) 

and an explicit relation for the attenuation constant in the vacuum: 

^- k i~0'- (40) 

An implicit solution for the surface wave frequencies is found by matching the solutions in 
Eqs.( l37"|) and f )38|) at z=0 using electromagnetic boundary conditions. The unique conditions 
are continuity of tangential H, E and continuity of normal D. When satisfied, the following 
dispersion relation results: 

(/3 + z47r X me ^)+e y /3 o = 0. (41) 

Note that the terms in Eq. fj4T|) do not have an odd multiple of wave vector k, and so 
the solutions should be reciprocal in the sense that u (k) = u(—k). We also see that the 
existence of surface modes strongly depends on the value of e y , since the solution of Eq.flJT]) 
for surface modes can only be found when the value of e y is negative. 

V. RESULTS FROM NUMERICAL CALCULATIONS 

We now illustrate the preceding theory for the material BaMnF4. Parameters appropri- 
ate for BaMnF4 were determined as follows. Measured values of the weak ferromagnetism 1 ^, 
M c = 18.347 Oe, and a canting angle 1 ^ of 3 mrad inserted into the relation M c = 2M S sin#, 
yield the magnetisation of the sub-lattices M s = 3.054 x 10 3 Oe. Following the approxima- 
tions by Holmes^ for the exchange field, He= 50T, and the relation He = XM S , the exchange 
constant is A=163.72. Using the measured value of the magnetic resonance frequency^, u r = 
3 cm -1 , in the relation u r = rfKiK — 2A)M^, we obtain the anisotropy constant ZT=0.337. 

The linear magnetoelectric coupling a is obtained by using the calculated spontaneous 
polarisation^ p o = 3.45 x 10 4 statC/cm 2 in Eq.©, yielding a = 1.42 x 10" 5 cm 2 /statC. 
The Ginzburg-Landau constant £1 and ( 2 are approximated by solving simultaneously Eq. ([8]) 
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(a)Dispersion relation. 
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(b)'Window' where surface modes exist. 
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(c)Narrow gap at Wj. (d)Narrow gap with bigger coupling. 

FIG. 2. Dispersion relation without external field. In (a) he dispersion relation without the external 
field is shown. The surface modes are indicated by"SP". The shaded regions represent bulk bands, 
which are limited by frequencies ui' m , co ez , uj oz , oj' ey and Wj. In (b) The "window" where the surface 
modes exist is shown. In (c) a narrow gap around cjj is expanded. In (d) the narrow gap is wider 
when the ME coupling is increased by a factor of 10. 



and (Ci + 3( 2 Po ) f = t ne soft phonon frequency along the spontaneous polarisation. 
With the inverse phonon mass 4 as / = 2.7 x 10 24 statC 2 /g cm 3 and transversal phonon 
frequency!- 9 u ey = 7.73 THz, this gives Ci=-10.528 g cm 3 /statC 2 s 2 and C2 = 1-934 x 10" 8 erg 
cm 5 /statC 4 . The suceptibility xt 1S calculated by using the transverse phonon frequency for 
z polarization 4 -, u ez = 33.7 cm -1 . 

The parameters C c and C a f in Eq.f l24"j) and f )33|) are converted to frequency units in cm -1 . 
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FIG. 3. Attenuation constant as a function of frequency. In (a) the real part of attenuation constant 
is shown for two values of wavevector in the absence of any external fields. The solid line represents 
k = cm -1 , while the dashed line corresponds to k = 250 cm" 1 . In (b) the imaginary part of 
attenuation constant is shown for two values of ME coupling. The solid line is for a = 1.42 x 10~ 5 
cm 2 /statC, and the dashed line represents the coupling a = 1.42 x 10 -4 cm 2 /statC. 



These become: 



C, 



Cat 



2 a 2 /i M 2 / 

71 (27TC) 2 

1 a/i MJ 



7T (27TC) 2 

with M s in A/m, a in m 2 /C, / in C 2 /kg m 3 and c in cm/s. 



(42) 



Solutions of Eq.f l36j) and ( l4"T|) are plotted in Fig. (2 (a) ) for the case with no applied fields 



present. There are two gaps in the bulk region which are created by the poles and zeros of 



Eq.( )36|) . First, as illustrated in Fig. (2(c)), a very narrow gap is located at the frequency U{ 
around 41 cm -1 , created by the magneto-electric interaction 4 -. This gap is associated with 
zeros in /(/i, e) = fi x e y — (47rx me ) 2 . The pole in the bulk modes is due to a zero of the dielectric 
constant e y . This gap is strongly dependent on the ME susceptibility, x me , and disappears 
when x me =0. The width of the gap is approximately proportional to {x me ) 2 - Thus the gap 



becomes wider with larger ME coupling. This increase is illustrated in Fig.(2(d)) where the 
ME coupling has been increased by a factor of ten. 

A second gap exists near the magnetic frequency w m =3 cm" 1 , and occurs at zero of 
f(fi) = fi x e y — (47rx me ) 2 at the magnetic frequency u) m . The other boundaries for the bulk 
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regions are determined by the attenuation constant. 

Using the attenuation constants, we obtain a narrow window between transversal and 
longitudinal phonon frequencies, u' ey and co^y, associated with the pole and zero value of 



e y (see Fig, 2(b)). In this figure, since the ME coupling is weak, hence the frequency ojl v 
is very slightly below the induced frequency Wj. Since the value of e y between these two 
frequencies is negative, surface modes can be obtained inside this narrow window. The 
surface modes start from the crossing between the lightline (u = ck) and the resonance 
frequency u ey and terminate at the longitudinal phonon frequency cu^y. The frequency ul v 
can be approximated as 



U L y = - < {U 2 ey + UJ 2 m + /) + (U 2 ey + UJ 2 m + /) - 4 {uj 2 ey UJ 2 m + UJ 2 m f) 



1/2 



(43) 



and is indicated in Fig, 2(d) Since the surface modes terminate at the longitudinal phonon 
frequency, the gap around uji does not influence the surface modes. 

Interestingly, from the expression for the attenuation constant of Eq. (I39[) . and the dis- 
persion relation of surface modes in Eq. (141 j) . it was only possible to satisfy the boundary 
conditions with a complex /3. The resulting mode is not a true surface mode but is instead 
a pseudo-surface wave^i. Because the solution for the attenuation constant in Eq. fl39l) is 
complex, the attenuation constant for the material sample consists of real and imaginary 



part as illustrated in Fig, 3(a) and 3(b) 
The imaginary part of /3 is 



, to 



-i47c X me -. 



(44) 



The positive value of the real part defines regions where the surface modes can exist. The 
existence of an imaginary part indicates that the solution in Eq.( l37|) is a psuedosurface 
mode, and not purely localized to the surface. Instead, energy "leaks" into the bulk. The 
wave is comprised of a localised component, that travels along the surface and decays into 
the material according to the real part of 0, and a component that travels into the material 
with wave number equal to the imaginary part of /?. 



Values for the imaginary parts of (3 are plotted as a function of frequency in Fig 3(b) 



for the case of no applied fields. Imaginary depends linearly on the magnetoelectric 
susceptibility and becomes large near the electric resonance frequency u ey . One can also 



see from Fig, 3(b), that the coupling directly influences the magnitude of "leakage". If the 
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coupling is large, then the ME susceptibility will also be large and thereby increase the 
imaginary part of j3. 

In the case where k » -, the attenuation constants for the material and vacuum regions 
can be approximated by e z (3 2 ~ e y k y and (3 Q ~ k y . The dispersion relation Eq.(l36l) then 
reduces to: 



e z e y 



1. 



(45) 



The surface modes require the permitivity e y to be negative, and the longitudinal phonon 
frequency polarized along z direction, u oz , is lower than u' ey . This means that the value of 
e z is positive in regions where surface modes exist, and the requirement in Eq. (1451) is never 
satisfied. Therefore, surface modes do not exist in the limit k » —. 

We now study the influence of external fields on the band structures. Results are shown 



in Figs 4(b) and 4(d) Results for an electric field value of 5x 10 V/m along the direction 



of spontaneous polarisation (but with zero magnetic field) are shown in Fig 4(b) The effect 



of the electric field is to shift the window where surface modes exist to higher frequencies. 
The upward frequency shift is due to the increase of spontaneous polarisation, which directly 
increases the phonon frequencies uj' and uj^ y . However, the change in uj Ly is smaller than 
the change in u' ey (which is due to the third term in Eq.( l43l ) and so the surface mode window 



is narrowed. The electric field increases the canting angle slightly, as shown in Fig, 4 (a) , and 
the effect on magnetic resonance is negligible. 



Results for a magnetic field of 10 T (with zero electric field) are shown in Fig, 4(d) The 



magnetic field increases the canting angle (see Fig, 4(c)) and shifts oo m to a lower frequency 



as shown in Fig 4(d) ) but the effects on the bulk bands are negligible. This shift can be 



understood if we consider the case where the ME coupling is neglected. Then the frequency 
Ume — > 0, and also Q me — > 0. In this case, the magnetic resonance frequency will take 
the form = O^fm + — w a(w« — 2u e x) cos 2 9. It can be seen from this expression 
that a magnetic field reduces the magnetic resonance frequency. We note that this effect is 
mentioned in Ref.— . 

The surface modes can be modified if the magnetic resonance frequency can be shifted 



into a surface wave "window", as shown in Fig, 5(b) and 6(b) In BaMnF4 where the electric 
and magnetic resonance frequency are well separated around 37cm -1 , it is very difficult to 
arrange the magnetic resonance frequency inside this window. This may be possible for a 
suitably prepared material (or artificially constructed material) whose frequency separation 
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FIG. 4. Influence external field in the canting angle and dispersion relation In (a)the external 
electric field E increase the canting angle with small value. In (b) dispersion relation is shown with 
an external electric field along the spontaneous polarisation E = 5 x 10 8 V/m. In (c)the external 
electric field H increase the canting angle. In (d), magnetic resonance frequency shift down when 
the external magnetic field, H Q = 10T is applied. The dotted line is the magnetic frequency in the 
absence of magnetic field. 



between electric and magnetic resonance is smaller. 



VI. PARAMETER EFFECTS ON SURFACE MODE PROPERTIES 

Lastly, we identify the key parameters affecting surface mode frequencies. In the first 
case, changing the phonon mass to / = 9 x 10 21 statA 2 s 2 /gcm 3 moves the magnetic resonance 
frequency to 1 cm" 1 above the the electric resonance frequency u ey . The result on surface 
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(a)Dispersion relation with close frequency of (b)The magnetic resonance inside the window 
magnetic and electric 

FIG. 5. Dispersion relation for a material with electric and magnetic resonances near one another 
in frequency. In (a) the dispersion relation without external magnetic field is shown. In (b), the 
magnetic resonance frequency is shifted into a surface mode window with the application of a large 
magnetic field, H = 12T. 



and bulk polariton bands is shown in Fig, 5 (a) The dielectric constant background has also 
been reduced to = 2.6, in order to widen the surface mode window. Application of 
an external magnetic field lowers the magnetic resonance frequency. Application of a large 
external magnetic of 12 T places the magnetic resonance inside the window as illustrated in 



Fig 5(b) 



Inside the window, the magnetic resonance splits the surface mode into low and high 
frequency branches for each direction of propagation. The properties of the upper part are 
similar to that discussed in the previous section. However, the lower branch terminates at 



the magnetic resonance frequency as illustrated in Fig. 5(b) In this case, the requirement 
that the dielectric constant e y should be negative for surface modes prevents both the upper 
and lower branches to exist in the region where k » — . 

In a second example, we consider if the exchange constant is -4000 and the anisotropy 
constant is 4. The ME coupling is also changed to 40 m 2 /C, which keeps the canting angle 



small. The dispersion relation at temperature 150 K is presented in Fig, 6 (a) As temperature 
is increased to 250 K, the polarisation will decrease while the magnetisation does not change 
significantly. Hence, the electric resonance goes below the magnetic resonance frequency. 
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FIG. 6. Temperature dependence with modified exchange constant, (a) Dispersion relation at 
T=150 K. In (b), the magnetic resonance frequency shifts down at 250 K 



Results are shown in Fig. 6(b) As in the first example, the magnetic resonance frequency 
then exists inside the surface mode window and a similar splitting of the surface mode 
branches occurs. 



VII. CONCLUSIONS 



We have shown how linear magnetoelectric coupling influences surface and bulk TM 
polaritons modes in a canted multiferroic. In the first place, a narrow restrahl region forms 
in the bulk mode band at a frequency near the longitudinal phonon frequency (along the 
spontaneous polarization). Reciprocal surface excitations can exist in this region. In general 
the surface exciations in this polarization are actually psuedosurface waves that are only 
partially localized to the surface. The imaginary part of the decay constant is proportional 
to the magneto-electric coupling. 

We have explored also possibilities of modifying the polariton band structure through 
changes in key material parameters. Such effects may possibly be realized in appropriate 
compounds, or in artifically layered heterostructures. If the magnetoelectric material has 
a magnetic resonance frequency above the electric resonance frequency, and the difference 
between those frequencies is not great, then it is possible to create surface excitations that 
are sensitive to temperature, electric fields, and magnetic fields. 
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